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The connectivity and tunability of superconducting quantum devices provide a rich platform to build quantum simulators and study 
novel many-body physics. Here we study quantum phase transition in a detuned multi-connected Jaynes-Cummings lattice, which 
can be constructed with superconducting circuits. This model is composed of alternatively connected qubits and cavity modes. 
Using a numerical method, we show that by varying the detuning between the qubits and the cavities, a phase transition from 
the superfluid phase to the Mott insulator phase occurs at commensurate fillings in a one-dimensional array. We study the phase 
transition in lattices with symmetric and asymmetric couplings, respectively. 
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1 Introduction 

Quantum simulation has become a frontier in quantum in¬ 
formation science after it was first studied in the pioneer 
work of Feynman m. State-of-the-art technologies in build¬ 
ing atomic and solid-state devices in the quantum limit en¬ 
able us to build quantum simulators with rich controllability 
and novel physics effects Among such efforts, super¬ 
conducting quantum simulation is being intensively explored 
in recent theoretical and experimental works. With super¬ 
conducting circuits that contain only one or two qubits, peo¬ 
ple have demonstrated Anderson localization, switch of the 
Chern number in a topological phase transition, and the si¬ 
multaneous coimling of a superconducting resonator to mul¬ 
tiple qubits |I3]Q. These experiments are enabled by tech¬ 
nological developments in controlling and coupling super¬ 
conducting devices with long decoherence times and high 
gate fidelity in the past decade GHU- Meanwhile, the 
experimental progress is paralleled by theoretical propos¬ 
als on building both digital and analog quantum simulators 
with superconducting circuits. The many-body models be- 
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ing studied include quantum spin systems lfT3UT8l . mod¬ 
els bearing novel topologic al prope rties lfT9l422ll . models of 
electron-phonon interaction II23II24I . and systems in high en¬ 
ergy physics Il25l427ll . One particularly interesting model 
is the so-called coupled cavity array (CCA) model that is 
made of arrays of optical or microwave cavities, in which 
each cavity couples to a nonlinear medium. It was shown 
that the CCA model can demonstrate Mott insulator (Ml)-to - 
superfluid (SF) phase transition for cavity pola ritons Il28l4^ . 
in analogy to the Bose-Hubbard model 04l[T5l . One ad¬ 
vantage of analog quantum simulators over general-purposed 
quantum computers is that they put less stringent require¬ 
ments on the quantum logic operations and are more robust 
against decoherence and leakage errors. 

In a recent work, we proposed a multi-connected Jaynes- 
Cummings (JC) lattice model, which consists of arrays of 
alternatively coupled qubits and cavity modes ll^ . This 
model can be constructed with superconducting systems, e.g., 
by connecting X-mon qubits and coplanar waveguide res¬ 
onators i). Using the exact diagonalization method, we stud¬ 
ied the quantum phase transition of this model in a finite- 
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size one-dimensional (ID) lattice. This system demonstrates 
phase transition from the incompressible MI phase to the gap¬ 
less SF phase by varying the ratio between the left and the 
right qubit-cavity couplings. Moreover, the system shows 
the novel feature of reentrance to the MI phase from the 
SF phase when the coupling ratio is further increased. This 
reentrant behavior originates from the left-right symmetry be¬ 
tween the couplings, and it distinguishes the multi-connected 
lattice model from the CCA. Note that effective cavity mode 
coupling and quantum magnetism were studied in intercon¬ 
nected qubit-cavity arrays that bear uniform or opposite cou¬ 
plings EHH). 

Here we study a detuned multi-connected JC lattice, ex¬ 
ploring the role of qubit-cavity detuning in the MI-to-SF 
phase transition. We show that the effective onsite nonlinear¬ 
ity strongly depends on the detuning. By varying the detuning 
from large negative to large positive values, the nonlinearity 
can be tuned from weak to strong, compared with the mag¬ 
nitudes of the qubit-cavity couplings. We use the exact diag- 
onalization method to find the many-body ground state of a 
ID lattice. Then, the single-particle density matrix at com¬ 
mensurate fillings is calculated. Note that even though ID 
bosonic systems do not possess long range order, the single¬ 
particle density matrix decays much faster in the MI phase 
than in the SF phase and still gives clear signature of the 
phase transition. We study the quantum phases for lattices 
with symmetric and asymmetric couplings, respectively. For 
symmetric couplings with the left and right couplings equal to 
each other, the phase transition demonstrates universal behav¬ 
ior that only depends on the ratio between the detuning and 
the coupling; whereas for asymmetric couplings, the phase 
transition shows more complicated behavior. 

The qubit-cavity detuning can be conveniently adjusted by 
applying a global dc magnetic field to the SQUID loop of 
the X-mon qubits HI. The magnetic field controls the energy 
level splitting of the qubits, and hence the detuning, over a 
very wide range. In comparison to varying the qubit-cavity 
couplings, a practical advantage of varying the detuning is to 
avoid additional circuit elements that could bring in serious 
complications in circuit design and operation. 

This paper is organized as follows. In Sec. |2] we present 
the Hamiltonian for a multi-connected JC lattice model in 
a ID array, and we analyze the effective onsite interaction 
and the effective hopping matrix element. The dependence of 
the onsite interaction on the qubit-cavity detuning is studied 
in detail. Sec. |3] starts with a description of our numerical 
method. We then present our results of the single-particle 
density matrix for both symmetric and asymmetric qubit- 
cavity couplings. Conclusions are given in Sec.|4] 


where crf’’^ denote the Pauli operators of the qubit, aj (a,) is 
the creation (annihilation) operator of the cavity mode, cJc is 
the angular frequency of the cavity, and the energy level 
splitting of the qubit. The qubit-cavity detuning A is defined 
as the frequency difference between the qubit and the cavity 
with A - a>c - cl>z- 

In addition, the qubit at site i couples to the cavity at site 
i - 1 with the interaction 

Ki^Si{aUo-- + ( 2 ) 

which is also in the form of the JC coupling. Combining the 
local Hamiltonian ([T]i and the interaction we can write 
the total Hamiltonian of the multi-connected JC lattice model 
as 

+ ( 3 ) 

i 

where the summation runs over the entire lattice. It is note¬ 
worthy that the total Hamiltonian is invariant under a left- 
right reflection of the lattice accompanied by an exchange 
of the couplings g[ and gr- This reflection symmetry is at 
the heart of the reentrant behavior studied in our previous 
work Multi-connected lattices with more complicated 
geometries in higher dimension can also be studied. For ex¬ 
ample, the qubits and cavities can be connected in a checker¬ 
board pattern in a two-dimensional lattice, as illustrated in 
Fig. mb). 


(a) 
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Qi-1 Rt-i Qi Ri Qi+i Ri+i 
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2 Multi-connected JC lattice 

2.1 Model 

The multi-connected JC lattice model is made of alternatively 
connected qubits and cavities 1^ . as illustrated in Fig.m^). 
In a ID configuration, each qubit couples to two neighboring 
cavities, one to the left and the other to the right hand side, 
with independent coupling constants gi and g^. We define a 
unit cell as one qubit and the cavity to its right hand side. 
The local Hamiltonian for the /-th unit cell has the form of a 
standard JC model with (Jt = 1) 

a)ca]ai + ^(r.+gr[a]a-J -Hcrta,), ( 1 ) 


Figure 1 Multi-connected JC lattice models, (a) A ID lattice made of 
alternatively connected qubits Qi and cavities /?/. Each qubit Qi couples 
simultaneously to Ri with coupling gr and to /?/_i with coupling gi. (b) 
A two-dimensional lattice with the qubits and the cavities connected in a 
checkerboard pattern. Here a unit cell consists of a qubit Q/j, a cavity R^., 
and another cavity R^^ The qubit Qij couples to four cavities: R^y R^y 
R^_^ y and with respective couplings. 

2.2 Nonlinearity and hopping 

The total Hamiltonian Ht possesses the key elements for a 
MI-to-SF phase transition: an effective onsite nonlinear inter¬ 
action and effective hopping matrix elements between neigh¬ 
boring sites. To illustrate this point, we study this model 
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in the strongly-asymmetric limit of gi gr (or similarly, 
gr ^ gi)- First, we consider the local Hamiltonian H'j^ at 
site i. The eigenstates of H'j^ include a zero-excitation state 
|0,) = |0, i) and the polariton doublets In,-, +) with 

\ni,+) = cos i) + sin y'”'“ 

3 3 

|u/,-> = sin Y|u/,i)-cos y|«, - 1,T), (4b) 

where | T. i) correspond to the spin up or down states of the 
Pauli operator (j~, the u, index on the right hand sides of the 
equations refers to cavity photon numbers, and the «, index 
in the eigenstates |n/, +) refers to the total number of exci¬ 
tations at site i that counts both spin flip and photon excita¬ 
tions. Here cos % = .^^(1 + sin ^ = .^^(1 - 

and -i- 4g^n,' is an effective Rabi frequency. The 

polariton doublets |n,, +,) are eigenstates with «,■ excitations. 
The eigenenergies can be written as 

e«„± = (h/ - 1/2)We + D.„J2 (5) 

and So, = respectively. The states |n/, -) have lower 

energy in the «, doublet, and are referred to as the lower po¬ 
lariton states. The energy difference Ae„; = - sq, is thus 

the lowest energy cost to produce «,■ excitations at site i. 

To analyze the effective nonlinearity, we use the follow¬ 
ing estimation. Assume that the energy of lower polariton 
states |«;, -) can be described by an effective Hamiltonian 

^fefF = ^pp]pi + {UI2)p]p]pipj, where Up is the single- 
polariton energy and U is the magnitude of the Hubbard on¬ 
site interaction. The energy cost to generate «, excitations is 
thus Ae„_ = iiiOJp + (f//2)«,(«; - 1), with the special case of 
Aei, = (jjp. With these relations, the Hubbard interaction can 
be written as U - (Ae„,^| - As,,. - Up)lni. Combining this re¬ 
sult with the eigenenergies given in Q, we can estimate the 
interaction U in terms of the detuning and the qubit-cavity 
coupling for selected n, value. Note that different from the 
Hubbard model, the effective U in our estimation depends on 
tij. The JC model cannot be exactly mapped to the Hubbard 
model, even though the nonlinearity is a convenient way to 
analyze the behavior of this system. For example, at zero 
detuning with A = 0, t/ = gr(l + Vul - Vn7+T)/n, for n, 
excitations, proportional to the coupling strength gr and de¬ 
creasing monotonically with the total excitation number u,. 
The dependence of U on «, is shown in the inset of Fig.|3 

In this system, the detuning can be an effective knob to 
tune the phase transition between the MI and the SF phases. 
Practically speaking, it is easier to control the qubit-cavity 
detuning than to adjust the qubit-cavity couplings in the lab¬ 
oratory. The detuning of superconducting qubits can be con¬ 
trolled with global external magnetic held over a very wide 
range HI. Below we study the dependence of the Hubbard 
U on the detuning. In Fig. we plot U as functions of 
the detuning for low-lying polariton states with n, = 1 to 
10 and a coupling strength of gr - 150 MHz. The nonlin¬ 
earity increases continuously with the detuning. In the limit 
of large detuning magnitude with |A| » go, we And that 
[/ ^ (A + |A|)/2n„ i.e., for large positive detuning, U A/n, 
and for large negative detuning, t/ —» 0. This indicates that 
for low-lying polariton excitations, the SF phase is more fa¬ 
vorable at negative detuning with a nearly vanishing effective 
interaction; and the MI phase is more favorable at large posi¬ 
tive detuning. 



A/g^ 

Figure 2 The effective nonlinearity U versus A/gr with gr = 150 MHz. 
Inset: U versus rij for A/gr = -2, 0 and 2, respectively. 


Assume the coupling gi between neighboring sites to be 
finite with gi <& gr- The coupling between adjacent 
sites can be treated as a perturbation to the eigenstates (l4al i 
and (Hbl) . Consider the matrix element of between state 
|A) = |n,-i, -)® |u/, -), which contains n/_i excitations at site 
i - 1 and n, excitations at site i, and state |.B) = |n/_i - 1, -) ® 
|u, + 1, -), which includes (u,_i - 1) excitations at site i - 1 
and (u, -H 1) excitations at site i. We And 


{B\crfai^i\A} = - sin Y cos x 

I- . — l 

Jui-i sm-sm- 

' 2 2 


I 7 bfl. j_i bn—1 

+ yrii-i - 1 cos —-— cos 


2 r 


(6) 


which can be viewed as a hopping matrix element for po¬ 
lariton excitations between neighboring sites, with a hopping 
amplitude proportional to the coupling gi. The total number 
of excitations in the lattice is preserved during the hopping 
process. Hence in the limit of very asymmetric couplings 
with gi « gr, the onsite coupling gr produces an effective 
Hubbard U and the weak nonlocal coupling gi produces a 
hopping term between adjacent lattice sites. Similar argu¬ 
ment can be applied to the opposite limit of g[ » gr, in which 
case gr becomes the hopping perturbation and gi yields the 
effective nonlinearity. 

Compared with the CCA, the multi-connected JC lattice 
model bears the distinct feature that the cavities do not inter¬ 
act directly, and as a result, the hopping of the excitations is 
not caused by direct coupling between cavity modes. In the 
CCA, local qubit-cavity coupling produces onsite nonlinear¬ 
ity and cavity coupling induces hopping between neighboring 
sites. Whereas in the multi-connected JC lattice, both onsite 
nonlinearity and hopping are caused by the qubit-cavity cou¬ 
plings. In particular, when the magnitudes of the left and the 
right couplings are comparable (gi ~ gr), the role of each 
coupling cannot be simply analyzed as either inducing the 
nonlinearity or generating the hopping. It is the interplay be¬ 
tween the left and the right couplings that give rise to the 
phase transition in this model. 
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3 Detuning induced phase transition 

3.1 Numerical method 

In this section, we study detuning-induced quantum phase 
transition in the multi-connected JC lattice model using the 
exact diagonalization method. For this purpose, we choose 
the basis vectors 

|i/r> = |ni,cri>®|n 2 ,cr 2 >---|«M,crM) (7) 

to construct the matrix representation of the total Hamilto¬ 
nian for a finite system with lattice size M and total excita¬ 
tion number N. Here n, is photon excitations, and cr, =X, "f 
is the qubit state at site i with J, (T) corresponding to 0 ( 1 ) 
spin excitation. The operator for the total excitations in the 
lattice is defined as iV = + crfcrj), including both 

photon and spin excitations. In our model, with the commu¬ 
tation relation [N, H,] - 0, the total excitation number is a 
good quantum number. Below we study the quantum phase 
of many-body ground state at commensurate fillings, i.e., the 
total excitation number is a multiple of the lattice size M. 
Hence, only these basis vectors that have N excitations are 
selected in the numerical calculation. Because the Hamilto¬ 
nian is a sparse matrix on the selected vector space, we use a 
Lanczos algorithm to solve the many-body ground state |G). 
Note that we have assumed the periodic boundary condition 
for the lattice to ensure translational symmetry. 

For a finite size system, the order parameter (G|a,|G), 
which is often used to characterize quantum phase transition 
in the thermodynamic limit, is always equal to zero because 
of particle number conservation. Instead, we use the single¬ 
particle density matrix defined as 

p,{i,j)^{G\a]aj\G)l{G\a]a,\G) ( 8 ) 

to characterize the phase transition in this model. From def¬ 
inition, is hermitian and semi-positive definite. Fur¬ 

thermore, due to the translational and reflectional symme¬ 
try of the multi-connected JC lattice under periodic boundary 
condition. Pi (i, 7 ) is also real, symmetric, and cyclic. Hence 
p\{i,i) only depends on |/ - 7 I, and hereafter we write the 
single-particle density matrix as p\{x) with x = I/- 7 I. For ID 
bosonic systems in the SF phase, pi(x) decays algebraically; 
whereas in the MI phase, pi(x) decreases exponentially to 
zero. Hence, pi(x) in the SF phase could have much larger 
value at finite x than that of the MI phase due to its slower 
decay. Note that due to the finite size effect, pi(x) in our cal¬ 
culation will not reach zero in the MI phase. It decreases to 
a small finite value that will be clearly distinguishable from 
that of the SF phase. As a result, we can verify the occurrence 
of the phase transition, even though we cannot accurately de¬ 
termine the position of the quantum critical point. 

In our calculation, we choose a lattice size M - % and a to¬ 
tal excitation number A = 8 , with the qubit-cavity couplings 
in a range of 0 - 300 MHz and the qubit-cavity detuning in the 
range of A/go 6 [-3,3]. Here go - 150 MHz is introduced 
as a unit of the detuning, with its strength comparable to the 
qubit-cavity couplings. 

3.2 Phase transition with symmetric couphngs 

First, we study the detuning effect on the quantum many- 
body phase of a symmetric multi-connected JC lattice model, 
in which the left and the right qubit-cavity couplings are equal 
to each other with gi-gr- The many-body phase of this sys¬ 
tem thus only depends on the ratio between the detuning and 
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the coupling gr- In Fig.0 we plot the single particle density 
matrix pi(x = Xmax) as a function of the detuning A at the 
maximal lattice distance Xmax and at selected values of the 
coupling gr- Note that for the lattice size M = 8 , Xmax = 4 
under the periodic boundary condition. 



Figure 3 The single-particle density matrix pi (x = Xmax) versus the unitless 
detuning A/go at selected couplings with gi = gr and go = 150 MHz. 

The numerical result shows that pi(Xmax) decreases mono- 
tonically with the detuning. For large negative detuning with 
-A » gr, Pi(xmax) ^ 1. At A = 0, Pi(xniax) ~ 0.75, which is 
still a large finite value of the order of unity. For large positive 
detuning with A » gr, however, pi(Xmax) —> 0. The system 
hence exhibits a transition from the SF to the MI phases as A 
increases from negative to positive values. This result shows 
excellent agreement with the nonlinearity given in Fig. |2] in 
which the effective nonlinearity increases from nearly zero to 
large finite value as the detuning increases. The phase transi¬ 
tion can be explained with an intuitive physical picture. For 
large negative detuning, ojc u^,, and the lower polariton 
states can be approximated as photon states with equal en¬ 
ergy spacing, i.e., having weak nonlinearity. The many-body 
ground state in the presence of finite qubit-cavity coupling 
is then a superfluid of the photons. For large positive detun¬ 
ing, on the other hand, Uc » and the lower polariton 
states are approximately the qubit states that are localized at 
its own site. The many-body ground state then resembles the 
MI phase where hopping is prohibited by large nonlinearity. 

Note that for symmetric couplings, the behavior of this 
system depends on the ratio /S.lgr only. Hence the curves for 
different gr values in Fig. |3] can be rescaled to become one 
single curve. In other words, the width of the transition re¬ 
gion between the SF and the MI phases scales with the cou¬ 
pling constant gr, with wider detuning region for greater gr 
value. It is also worth noting that pi(Xniax) at A = 0 is inde¬ 
pendent of the coupling gr with the system in the SF phase 
regardless of the coupling strength. 

3.3 Phase transition with asymmetric couplings 

Now we consider the detuning effect on the phase transi¬ 
tion of asymmetric systems with distinct qubit-cavity cou¬ 
plings g! and gr- We calculate pi(Xmax) as a function of A 
for selected pairs of couplings as shown in Fig. |4] The cou¬ 
plings are chosen in the experimentally approachable regime 
of gir e [0, 300] MHz and satisfy gi + gr ^ 300 MHz. The 
numerical result shows the generic feature that the single¬ 
particle density matrixpi(Xniax) decreases with the increase of 
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the det unin g. In contrast to the case of symmetric couplings 
in Sec. 13.21 the ratio of the qubit-cavity coupling gilgr plays 
an important role in the many-body phase of the ground state. 
For instance, at {gi,gr) = (5,295) {gilg,- = 0.017),pi(vniax) ~ 
0 in the entire region of A/.go 6 [-3, 3], which indicates that 
the system is in the MI phase. Whereas at {gi, gr) — (50,250) 
igilgr = 0.2), pi(vniax) > 0.8 for ^|gQ = -3, dearly indi¬ 
cating that the system is in the SF phase. Hence depending 
on the ratio between the couplings, the phase transition can 
occur at either negative or positive detuning. 





Figure 4 The single-paiticle density matrix p\ (x = Xmax) versus the unitless 
detuning A/go at selected pairs of couplings and go = i50 MHz. The cou¬ 
plings are (g/,g,) = (0,300), (5,295), (25,275), (50,250), (100,200), and 
USO, 150) MHz, respectively. 


4 Conclusions 

To conclude, we studied the quantum phase transition in the 
many-body ground state of a ID multi-connected JC lattice 
model. Using a numerical method, we showed that by vary¬ 
ing the qubit-cavity detuning, the SF-to-MI phase transition 
can be observed in this system. The quantum critical point 
not only depends on the magnitude of the couplings, but also 
depends on the ratio between the left and the right couplings 
of the lattice. Our results show that the multi-connected JC 
lattice model could be a rich system to explore the many-body 
physics of cavity polaritons. 


This work is supported by the National Science Foundation under Award 
Number 0956064. 
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